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correlations between two disjoint blocks of a quantum critical system, to show that the 
structure of the causal cones characteristic of MERA, requires a transition between two 
different regimes attainable by changing the ratio between the size and the separation of 
the two disjoint blocks. We argue that this transition in the MERA causal developments 
of the blocks may be easily accounted by an AdS f j + 2 black hole geometry when the mutual 
information is computed using the Ryu-Takayanagi formula. As an explicit example, we 
use a BTZ AdS3 black hole to compute the MI and the quantum correlations between two 
disjoint intervals of a one dimensional boundary critical system. Our results for this low 
dimensional system not only show the existence of a phase transition emerging when the 
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transition when finite temperature and finite size effects are taken into account. 
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1 Introduction 

Entanglement entropy (EE) is by now regarded as a valuable tool to witness the amount 
of entanglement in quantum field theories and many body systems. By partitioning a 
given system S into two complementary sets A and A such that S = A U A, the reduced 
density matrix p A (i.e, the density matrix for an observer accessing only the degrees of 
freedom of the subsystem A), is obtained by tracing the full density matrix p over the 
degrees of freedom contained in A i.e, p A = Tr The EE accounts for the amount of 
quantum correlations between the complementary regions A and A and is defined as the 
von Neumann entropy of p A , 

S A = -Tr(p A log p A ) . (1.1) 

A standard approach to compute the entanglement entropy makes use of the replica 
trick [1-3]. The replica trick may be applied when the density matrix for the full system 
is represented by a path integral (as in the vacuum or in a thermal state); then, one can 
rather easily obtain the EE (1.1) of the subsystem A from the knowledge of, 

S A = -^Tvp A \ n=1 . (1.2) 

In [3], it has been shown that, for d = 1 quantum critical models, Tr//^ = c n (£ A /e)~^ c ^ G ^ n ~ 1 ^ n \ 
where i A is the length of the interval A, c is the central charge of the conformal field theory 
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(CFT) describing the given system at criticality and e is an ultraviolet cutoff. Using (1.2), 
one obtains that the EE is given by, 



S *=fa(j) + Sl , (1.3) 

where s\ is a non universal constant. 

Using an alternative approach based on holography Ryu and Takayanagi (RT) derived a 
celebrated formula yielding the EE of the region A provided that the (boundary) conformal 
field theory describing the critical system admits an holographic gravity dual [4, 5]. In the 
RT approach, the EE is obtained from the computation of a minimal surface in the dual 
higher dimensional gravitational geometry (bulk theory); as a result, the entanglement 
entropy £4 in a CFTd+i is given by the celebrated area law relation, 

Area( 7 ^) 

N 

where d is the spatial dimension of the boundary CFT, 7^4 is the d-dimensional static 
minimal surface in AdS^+2 whose boundary and area are given by dA and Area(7^i), re- 
spectively. G^ +2 ^ is the d + 2 dimensional Newton constant. The RT proposal is physically 
appealing since looking for the minimal surface 7,4 separating the degrees of freedom con- 
tained in region A from those contained in A amounts to search for the severest entropy 
bound on the information hidden in the AdS^+2 region related with A. For d = 1, eq. (1.4) 
becomes [4], 

Length^) 

Although the RT formula has not been rigorously proven its validity is supported by 
very comforting evidence. 1 For instance, one may show [8] that it provides a simple tool 
to prove the strong subadditivity of EE, i.e given two regions A and B, 

Sa + Sb > Saub + Sahb ; (1.6) 

furthermore, eq. (1.4) together with (1.6) may be used also (at least in the context of 
strongly coupled gauge theories, i.e at a t'Hooft coupling A > 1) to derive the concavity 
property of coplanar Wilson loops defined on curves Ca = dA and Cb = dB lying in the 
same two dimensional plane. Namely, 

(W(C A )) (W(C b )) < (W(C A ub)) (W(C AnB )) , (1.7) 

where Caub = d(A U B) and Cadb = d(A n B). To derive (1.7) one only needs to note 
that, from the Maldacena conjecture [9], the expectation value of a Wilson loop defined 
along a curve C is related to the area of the minimal surface 7 bounded by C by, 

(W{C)) ~ exp(-\/AArea(7)) . (1.8) 



1 See [6, 7] for some interesting attempts to derive it. 
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If in (1.8) one takes C = Ca = dA and 7 = 74, using (1.4) one can establish, up to 
a constant that ~ — log (W(dA)). Similar arguments yield Sb ~ — log (W(dB)), 

Saub log (W(d(A U B))) and S A nB log (W(d(A n £))). As a result, using (1.6), 

one gets (1.7). 

The minimal curves used in the RT formula, allow also to compute the two point 
functions of conformal primary operators of CFT^+i with an holographic gravity dual that 
is an asymptotically AdS^+2 space-time. The holographic computation of the correlation 
functions of these operators yields [10], 

(P(xi)0{xj)) ~ exp(-ALength( 7ii )) , (1.9) 

where A is the operator scaling dimension and 7^ is minimal curve in the bulk geometry 
connecting the boundary points Xi and Xj. 

Very interesting issues [11], [12] arise if one regards A as the union of several disjoint 
regions A = D{Ai and A as its complement. In the simplest case one may consider two 
disjoint blocks A and B such that A = A U B. In the analysis of those situations it is most 
convenient to compute the mutual information (MI) between regions A and B, which is 
defined by 

I(a-.b) = $a + Sb - Saub ■ (1-10) 

MI measures the amount of correlation (classical and quantum) between the spatially 
disconnected regions A and B and acts as an upper bound on the quantum correlations 
between operators defined in those regions [13]; namely, 



((OaOb)-(Oa)(Ob))' 



The correlators (OaOb) as well as I(a-.b) f° r two spatially disconnected regions disclose 
relevant information about the spatial distribution of entanglement in a given state of the 
system. However, for two disjoint blocks, neither the MI nor the quantum correlation 
functions happen to be a proper measure of the entanglement since A U B is not a pure 
state. A true measure of entanglement, requires the computation of negativity [14] which 
is a quite challenging task using field theory methods. 2 

By means of the replica trick, the computation of I(a-.b) requires the knowledge of 
Saub = — Tt(/?aub log/>AuB) with pavjb = Tr^p, for which very little is known so far. 
For two spatially separated regions A and B, the only exact result for Saub, has been 
obtained for free massless fermions in two dimensions [17, 18] but it remains unknown in 
its general form for other physically relevant theories such as the free compactified boson 
[12]. In a recent paper [12], Trp AuB for two disjoint intervals A and B - of length / 
(\u\ — vi\ = \u2 — U2I = separated by a distance d = \v± — U2I- has been computed 
yielding, 

Tr, 5ufl = 4 ( l*-«.lh«--»l ) §< - 1/ " > Mx) . (1 . 12) 

\«1 - 1*2 fl - V 2 \\Ui - V2\\Vi - U2\J 



2 See, for instance, [15, 16] for a discussion of this issue and some numerical examples. 
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with x being the conformal four-point ratio defined as 



\ui - vi\\u 2 - v 2 \ I ,, 1Q . 

x = -. n r = — —7 . (1-13) 

\ui - U2\\vi - v 2 \ [I + dy 

The function J- n {x) depends explicitly on the full operator content of the theory and is, 
of course, model dependent. However, the analytic continuation of J- n (x) to n = 1 in eq. 
(1.12) is hard to attain and this makes the computation of I^a-.b) between disconnected 
regions a rather difficult task within this approach. 

In a recent work [19], using the RT formula for d = 1 quantum critical system, it 
has been predicted the occurrence of a phase transition probed by the computation of the 
MI between two disjoint intervals of the boundary CFT^+i; namely, as the conformal four 
point ratio crosses a critical value the MI vanishes. Using exact methods, the vanishing of 
the MI has been confirmed to occur also for the critical XX spin chain [20]. This result is 
quite surprising from a quantum information point of view since, when the MI vanish, the 
PAub factorizes into pab = PA®Pb, implying that the two blocks are completely decoupled 
from each other and, thus, also the entanglement should be rigorously zero. In [19] it has 
been pointed out that, 

f°> x < 1/2 , 

IiA:B)(x) = \/ , (1.14) 



(c/3) log (*/(!- z)) , s>l/2 



where x is the conformal four point ratio defined in (1.13). Equation (1.14) states that 
I(A-.B) = for x < 1/2 and it has a discontinuous first derivative at xq = 1/2. As argued 
in [19], the discontinuity in the first derivative of the MI occurs since the shape of the 
geodesies (i.e., of the minimal surfaces in the bulk connecting the two disjoint intervals 
of the boundary critical system) changes, as x varies, due to the switching between two 
saddle points of the Euclidean action [6, 19] much similar to the one observed in [21]. 

In this paper, inspired by the analysis carried in [22-24], we exploit the holographic 
structure of the Multiscale Entanglement Renormalization Ansatz (MERA) tensor net- 
works, to analyze the correlations between disjoint blocks of a critical system described by 
a (d + 1) dimensional conformal field theory lying at the boundary of an asymptotically 
AdSd+2 spacetime. In order to get an hint on the pertinent ansatz for the metric to be used, 
we observe that, when computing the quantum correlations between two disjoint blocks 
of a boundary quantum critical system, the structure of the causal cones characteristic of 
MERA [25, 26] implies the existence of two different regimes attainable by changing a pa- 
rameter depending on the ratio between the size and the separation of the disjoint blocks. 
We argue that this transition may be accounted by an AdS^ + 2 black hole geometry and 
use the RT formula to compute the MI between the two disjoint regions of the boundary 
critical system. As an explicit example, we use a BTZ AdS3 black hole to compute the MI 
and the quantum correlations between two disjoint intervals of a one dimensional bound- 
ary critical system: here, our analysis not only confirms the existence of a phase transition 
emerging when the conformal four point ratio reaches a critical value but also provides a 
rather intuitive entropic argument accounting for the source of this instability. Finally, we 
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investigate how the holographic computation of the MI between two disjoint blocks may 
be affected by finite size (and temperature) effects. Besides its appealing beauty, we feel 
that a remarkable merit of the holographic approach is that it can help in establishing 
fruitful connections between the phase transition analyzed in [6, 19] and analogous phase 
transitions exhibited by disconnected operators such as the one occurring for disconnected 
Wilson loops found in [27]. 

The paper is organized as follows: in Section 2, we review the MERA induced AdS / CFT 
duality [22, 24] and analyze its relationship with the RT holographic formula [4, 5]; there, 
we argue that, when considering two disjoint blocks of the boundary CFT describing the 
critical system, the MERA induced AdS/CFT duality leads rather naturally to the emer- 
gence of an AdS black hole as the relevant space time metric in the dual bulk space. In 
Section 3 we briefly review the geometrical properties arising when the space time metric 
in the bulk is described by an AdS3 BTZ black hole; there we point out also how a BTZ 
black hole metric in the MERA induced dual AdS3 space easily accounts for the finite tem- 
perature corrections to the EE. In Section 4 we use the RT formula [4, 5] to compute the 
MI and the quantum correlations between two disjoint intervals in the CFT2 dual to the 
AdS3 BTZ geometry; there we show that the RT formula, when computed using the AdS3 
BTZ geometry, naturally accounts for the phase transition discovered in [19] and provide 
an entropic argument accounting for the emergence of this instability. Finally, in Section 5 
we summarize our results. In the appendix A we use our approach to compute the MI and 
quantum correlations between disjoint intervals of the boundary quantum critical system 
when the metric of the MERA induced AdS3 space is described by a spinning BTZ black 
hole. 

2 MERA induced AdS/CFT duality 

In [22] , it was firstly observed that MERA [25] may give rise to a realization of the AdS / CFT 
correspondence [9]. This observation has been subsequently developed in [24]. MERA is a 
real space renormalization group technique based on a series of consecutive coarse-graining 
transformations reducing the amount of entanglement in a block of lattice sites of a critical 
system before truncating its Hilbert space. Namely, by renormalizing the amount of entan- 
glement in a given system, the MERA procedure controls the growth of the sites Hilbert 
space dimension along successive scaling transformations. This entanglement renormaliza- 
tion procedure may be encoded in a tensor network arranged in a set of different levels 
{tffclj^Q accounting for the consecutive renormalization steps and, for quantum systems 
at criticality, it shows a characteristic fractal structure. The tensor network implements a 
renormalization group transformation which is local in space and scales local operators into 
local operators. Furthermore, using MERA, it has been shown that quantum correlations 
in the ground state of one and two dimensional critical quantum many body systems, could 
be arranged in layers corresponding to different length scales i.e to different steps in the 
renormalization process. 

As pointed out in [22-24], the entanglement structure in a quantum critical many body 
system, defines a higher dimensional geometry via the renormalization process described 
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by the scale invariant MERA tensor network. The emerging geometry can be engineered as 
follows: all the sites in the MERA tensor network are arranged in layers, each representing a 
different scale (coarse graining renormalization step). As a result, besides the coordinates 
labelling the position and the time t, in MERA, one may add a "radial" coordinate z 
labelling the hierarchy of scales. Then, the higher dimensional geometry defined by MERA 
may be usefully visualized by locating cells around all the sites of the tensor network 
representing the quantum state: these cells are unit cells filling up the emerging "bulk" 
geometry and the size of each cell is defined to be proportional to the entanglement entropy 
of the site in the cell. As a result of this procedure a gravity dual picture of the bulk emerges 
quite naturally from the entanglement of the degrees of freedom of the critical system lying 
on the boundary [28]. 

The discrete geometry emerging at the critical point is a discrete version of Anti de 
Sitter space (AdS) [22, 24]. For a one-dimensional quantum critical system with a space 
coordinate labelled by X, the continuous isometry w — > w + a, X — > e a X of the metric 

ds 2 ~ dw 2 + e~ 2w dX 2 , (2.1) 

is replaced by the MERA's discretized version, w — > w + k, X -> 2 k X or X -> 3 k X 
depending on the binary or ternary implementation of the renormalization algorithm [26] . 
The analog of w = log z in the tensor network is simply the variable labelling the number 
of renormalization steps carried out by the MERA algorithm. 

When considering a continuous version of MERA [29] , the discrete AdS-like geometry 
given by (2.1), approaches its continuous version i.e the 3-dimensional AdS space with the 
scale invariant metric, 

f 2 

ds 2 = (-dt 2 + dz 2 + dX 2 ) . (2.2) 

In (2.2), i is a constant called the AdS radius; it has the dimension of a length and it is 
related with the curvature of the AdS space. With this choice of the space time coordinates 
the one dimensional quantum critical system lies at the boundary (z = 0) of the bulk 
geometry. 

2.1 Ryu-Takayanagi formula in the MERA induced AdS/CFT correspon- 
dence. 

There is a striking relationship between the RT formula and the computation of the en- 
tanglement entropy in MERA. Using MERA, the reduced density matrices and hence the 
quantum correlations are determined by the structure of the causal cones [25]. The causal 
cone CC(B) of a block B of I sites of the boundary critical system, is determined by grouping 
- following all the levels of the MERA tensor network- all the renormalizing operators and 
sites which may affect the sites in the block S. As a result, to compute the entropy Sq of 
the block B it is necessary to trace out any site in the bulk geometry defined by the tensor 
network which does not lie in the CC{B) of the block. 

The boundary of the CC{B) is a curve 75 in the MERA induced AdS higher dimensional 
geometry. The length of 75 is, by definition, the sum of the entropies of all the traced out 
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sites [22, 25], and thus provides an upper bound for the entropy Sq of B [24], 

S B < Length( 7B ) . (2.3) 

The close relationship with the geometrical RT formula comes about when one realizes 
that 7b can be regarded as the minimal curve of RT, since it counts the minimal number of 
sites which must be traced out in the coarse graining process. Indeed, the minimal curve 
7B in an optimized scale invariant MERA network [24] has proven to saturate the bound 
given in (2.3); this has been confirmed by explicit computation in one dimensional critical 
systems, where it has been shown that Sb ~ § log£ [25]. Since 70 arises as a boundary of 
the CC(B), it can be interpreted as an holographic screen which optimally separates the 
region in the bulk described by the degrees of freedom of B, from its complementary region 
B. 

2.2 Quantum correlations between disjoint blocks from MERA. 

In order to use MERA for computing the two point correlation functions, one should first 
observe that, the CCs of two operators located at points s\ and S2 of the boundary critical 
system always grow (i.e the number of sites inside a CC(sj) at MERA level Wf~, is always 
bigger than the number of sites at level Wk-i), since k increases as one gets deeper into the 
bulk geometry defined by the tensor network (2.2). As a result, there is a level u>* where 
CC(s\) and CC{s2) overlap. When the CCs overlap, the operators defined on the boundary 
are correlated with an algebraic decaying functional dependence [25]. At variance, the CCs 
of operators defined on finite size disjoint blocks of the boundary critical system, tend to 
exponentially shrink along the "coordinate" w labelling the MERA level [25]. 

As a result, for two disjoint blocks A and B of the same size I, two situations may 
occur (see Fig 1) depending only on the distance d between the blocks: 

i) after u>h ~ log I renormalization steps, the CC(^4) and the CC(B) shrink to one after 
they overlap (Fig 1 top). Here one expects that the correlations between the two blocks of 
the boundary critical system decay algebraically. 

ii) after vjh ~ log? renormalization steps, the CC{A) and the CC(B) shrink to one 
without overlapping (Fig 1 bottom). Here one should expect that the correlations decay 
exponentially with the "distance" between the two blocks. 

It is easy to convince oneself that, if one defines w* ~ logc?, is realized when 

wh > (wh < w*). 

In the following of this paper, we make the ansatz that an holographic dual spacetime 
that may efficiently account for these two distinct behaviours of the casual cones, is given 
by an AdS^+2 black hole geometry of radius zh = I {wh = log /) when the MI between two 
disjoint blocks A and B is computed by means of the RT formula. To support our ansatz 
we explicitly compute the MI and the quantum correlations between disjoint blocks of a 
one dimensional quantum critical system using the RT formula for the AdS3/CFT2 corre- 
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Figure 1. Top: Schematic representation of MERA CC for two disjoint finite intervals A and B 
when the separation between them allows for overlapping after ~ log d renormalization steps. 
The overlap occurs before the causal cones shrink to one (in our representation, when causal cones 
stabilize their width after wh = log I renormalization steps). Bottom: Schematic representation 
of MERA CC for two disjoint finite intervals A and B when the separation between them does 
not allow for overlapping after w* ~ log d renormalization steps. The curve 7 that goes through 
the links between the nodes of the MERA network is the minimal curve separating the CC(A) and 
CC(B) from the traced out sites in the MERA bulk geometry (2.2). 

spondence with the bulk metric given by a AdS3 BTZ black hole. Under these assumptions, 
from equation (1.5), the MI reads 

!{A-.B) = — [Length^) + Length^) - Length (jaub )) , (2.4) 

with ja, 1b and 7aub being geodesic curves in the BTZ black hole background [30]. As 
we shall see in the following sections, a computation of the MI carried using this approach 
supports- from a different point of view- the results obtained in [19]. 

3 The BTZ black hole 
3.1 BTZ black hole solution 

Bahados, Teitelboim, and Zanelli (BTZ) showed that (2+l)-dimensional gravity has a black 
hole solution, the BTZ black hole, differing from the Schwarzschild and Kerr solutions 
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mainly in that it is asymptotically anti-de Sitter rather than asymptotically flat. The 
BTZ solution is clearly a black hole: it has an event horizon and (when rotating) an inner 
horizon, and it exhibits thermodynamic properties much like those of a (3+l)-dimensional 
black hole [30]. 

The BTZ black hole may be obtained by orbifolding AdS3 through SL(2, C) identifica- 
tions [31] and is a solution of pure gravity in three dimensions with a negative cosmological 
constant described by the Einstein-Hilbert action supplemented by boundary terms [31], 

I = Jd 3 x^ (R - A) + / bndy . (3.1) 

In the following we use the Euclidean signature, and use the notation of Misner, Thorne, 
and Wheeler [32]; as a result, r = 1/z so that the boundary is now located at r — > oo. A 
simple solution of the equations of motion is just the AdS3 spacetime, 

ds 2 = (1 + r 2 /£ 2 )dt 2 + ^ + r 2 # 2 (32) 

1 + r z /£ z 

AdS3 has maximal symmetry, with the isometry group being SL(2, C) = SL(2, R)^ x 
SL(2, R)r. 

A more general one-parameter family of solutions is provided by the non-rotating BTZ 
black hole of mass M [31], 

ds 2 = —p^dt 2 + ^^dr 2 + r 2 dtf , (3.3) 

describing an AdS black hole with an event horizon located at r = r + = i\J 8GM at a tem- 
perature T = r + /2ir£ 2 ; of course, for large r, the solution (3.3) asymptotically approaches 
AdS 3 . 

The metric of a rotating BTZ black hole of mass M and angular momentum J is given, 
instead, by 

ds 2 = ^"^I 2 "^^ ^ + fjl _ 2 , dr° + r 2 (# + i^dtf . (3.4) 



2^2 ^^"2 y2 ^ ^y»2 ^«2 ^ ^^-a2 



with 



r± 




(3.6) 



Now, the event horizon is located at r = r + with r + > r_ and r_ being the inner Cauchy 
horizon. Rotating BTZ black holes have been recently shown to be relevant in investigations 
of helical Tomonaga-Luttinger liquids [33] . 
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3.2 Dual CFT to the BTZ solution 

The boundary of asymptotically A0IS3 spacetimes is a two dimensional torus on which one 
can define a dual CFT with its conformal symmetry being generated by two copies of the 
Virasoro algebra acting separately on the left and right moving sectors. As a result, the 
CFT splits into two independent sectors at thermal equilibrium with temperatures, 

r + + r_ 



Tr 



(3.7) 



2vr£ 2 ' " 2tt£ 2 ' 
The mass M and the angular momentum J in the rotating BTZ black hole geometry are 
related to the Virasoro charges of the dual CFT on the boundary by 



u 



c 

24 



1 



(M£ + J), L 



c 

24 



1 



lQG y ' 24 16G 

with c = c given by the Brown-Henneaux holographic relation [34] 

3f 



{Ml - J) . 



2G 



(3) • 

V 



In the non rotating BTZ metric (3.3) one has 



C 

24 



U 



c 

24 



16G£ 



(3.8) 



(3.9) 



(3.10) 



It is easy to prove [31] that, for the AdS3 metric, (3.2) Lq = Lq = 0; this is just a conse- 
quence of its invariance under the SL(2, R)l x SL{2, R)# group of isometries generated 
by L 0j ±i and L ,±i- 

3.3 Geodesies in the BTZ geometry 

The RT formula uses the spacelike geodesies in a given metric. For the BTZ black hole 
these geodesies are well known. The length of the geodesies connecting two points Xi and 
Xj separated by a distance \xi — Xj\ and located at the boundary of the AdS3 space whose 
metric is described by a BTZ Black Hole (3.3), can be written as [35], 



L(xi,Xj) = 2£log 



/3 



sinh 



TTC 



TT Xi 



ft 



(3.11) 



with ft = 2irl 2 /r + and e the regularizing boundary cut-off. 

Using the RT formula (1.5) and (3.9), one gets the well known formula [3] for the EE 
of a single connected block A of length £a = \x± — X2I from the BTZ geometry; namely, 
one has that 



Sa 



(c/3)log (£ A /e) 



L(x 1 ,x 2 ) 
4Gg) 



c. \ft . , {tt£a 
- log — sinh — 

3 7T6 \ p 



r + — > {ft — > 00) 



(3.12) 



{7Tc/3){£ A /ft) + f log(/3/27re) , r+ ^ 00 {ft -> 0) 



- 10 - 




Figure 2. Geodesic used in computing the entanglement entropy of a region of length I = \x\ — x%\ 
in the AdS Black Hole geometry. Left: 7 does not approach the horizon (dotted line). Right: 7 
wraps around the horizon. 



As expected, from eq. (3.12) one recovers the logarithmic dependence only in the zero 
temperature limit (i.e., when the size of the interval is small in comparison with the 
distance of the horizon from the boundary); indeed, in this limit, the BTZ geodesies stay 
close to the boundary and only probe the asymptotic form of the AdS3 BTZ geometry. 
(Fig. 2 Left). At variance, when the size of the simply connected region A is bigger than 
the distance of the horizon from the boundary, the BTZ geodesies probe the black hole 
horizon extending tangentially to it; this induces the linear correction to the EE which, for 
a single connected region A, describes now a thermal state at temperature T = 1//3 (Fig. 2 
Right). 

For a rotating black hole, the geodesies are given instead by 



L(xi, Xj) = 21 log 



PlPr 
7r 2 e 2 



sinh 



7T\Xi X j . 

1 smii 



IT X j X n 



(3.13) 



where Pl,r = 1/^l,r- As a result one gets that 



i log 



— sinh 



TTC 



+ |log 



— sinn 
ire 



PR 



(3.14) 



Equation (3.14) factorizes into left and right moving sectors as expected from the left-right 
decoupling of the CFT2. 



4 Holographic Computation of Quantum Correlations and Mutual In- 
formation for two disjoint intervals 

In this section we derive both the MI and the quantum correlations between two disjoint 
intervals of a one dimensional critical system described by a CFT2 located at the boundary 
of the AdS3 space. We assume in the following that the disjoint intervals A and B have 
equal size I and are separated by a distance d. Namely, we take A = [u±, v±] , B = [1*2,^2] 
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with \u\ — v\\ = \u2 — V2\ = I and \vi — U2\ = d (see Fig. 3). We shall see how, in both 
computations, one can find a critical value of a pertinent parameter at which there is a 
transition between two very different behaviors. 

For the two disjoint intervals A and B the holographic computation of the MI requires 
to determine the minimal curve in the bulk homologous to A U B. In the MERA induced 
AdS/CFT correspondence, the curve "Yaub is generated by tracing out the bulk sites lying 
outside the CC(A) and CC(B) and is an holographic screen for the entropy contained in 
AU B. As a result, for generating this holographic screen, there are- just as in [19]- only 
two possible options given by (Fig 3 Left) and 7^,b (Fig 3 Right), respectively: the 

curve Jaub (Taub) corresponds to the overlapping (non-overlapping) configuration of the 
causal cones CC(A) and CC(B) depicted in Fig 1. Namely, Jauw describes a situation in 
which the two intervals are enough separated so that Laub = Li(ui, v\) + £2(^2^2)1 while 
^AuB ' describes a situation where the separation between the intervals is so small that the 
minimal curve of the region A U B, connects the inner and outer boundaries of the two 
regions so that Laub = Li(ui,v 2 ) + L 2 (vi,U2). 

(dis) 

Of course, when Jaub * s use d m the holographic computation of the MI, the MI 
vanishes as a consequence of eqs. (1.5) and (2.4). At variance, when one uses JauB' the 
holographic computation of the MI (2.4), using as the metric of the AdS3 bulk space the one 
corresponding to a BTZ black hole with the horizon located &t z+ = I from the boundary 
i.e f3 = 2-KZ+, yields (3.9) 

sinh (ttT\ui — v±\) sinh (ttT\u2 — ^2!) 



!(A:B) = 3 log 



(4.1) 



sinh (ttT\ui — V2D sinh (irT\vi — U2I) 
with T = l/p. 

One sees that Iia-.b) m ec i- (4-1) equals zero when a certain ratio between I = \u\— v\\ = 
\ u 2 — V2\ and d = \v\ — U2\ is reached; namely, one sees that the MI, when computed us- 
ing the BTZ black hole as the metric of AdS3, vanishes at a value of the conformal four 
point ratio given by xq ~ 0.53. This is in agreement with the result of [19]. However, an 
advantage of the MERA induced AdS/CFT correspondence lies in the fact that one can 
provide a rather intuitive entropic argument accounting for the use of either one of the 
two minimal curves depicted in Fig. 3 when performing the holographic computation of MI. 
Indeed, since the length of the curves Javb an d Jaub are ~ °y definition- the sum of the 
entropies of all the traced out sites, the transition between the two behaviors of MI occurs 
when the separation between the two (equal size) disjoint blocks A and B is such that the 
entropy due to the the tracing out process yielding equals the entropy due to the 

tracing out process yielding Jaub- 

A similar transition is found also in the computation of the quantum correlations 
between two primary operators 0{xa) and 0{xb) (xa £ A and xb 6 B) with conformal 
dimension A defined in the CFT^ + 1 describing the boundary critical system. This should 
be expected in view of the bound (1.11). The AdS/CFT correspondence implies [9, 36] 
that 

{0{x A )0(x B )) ~ e - mL( - XA ' XB ^ , (4.2) 
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La u b = + L 2 





U, 



u, 



Figure 3. Minimal curves used in the holographic computation of Saub and I(a-.b) f° r to disjoint 
intervals A and B. 



where A ~ ml and L(xa,xb) is the length of the shortest geodesic connecting the boundary 
points xa and xb- Using for L(xa,xb) the expression given in (3.11), one easily gets, 




x A - x B 



(0(x A )0(x B )) 

-2A 



ttT 



2A 



2A 



exp (—2ttTA\xa — xb\ 



sinh (itT\xa — xb\ 
z + > \xa — xb\ 

, z + < \x A - x B \ 



(4.3) 



From (4.3) one sees that there is a change from an algebraic to an exponential decaying 
behavior of the two point quantum correlation function and that the transition between 
the two regimes occurs when (0(xa)0(xb)) ~ e~ A . When this happens, one has that 



2-7T T \ xa — xb\ 



\xa - xb\ 



1 



(4.4) 



which defines the value of the parameter fi = \xa — xb\/1 &t which this transition occurs 
[37]. 

Eqs. (4.1) and (4.3) are derived for infinite systems when the central charge c —> oo. 
However, one may be interested in the behavior of the MI and of the quantum correlations 
in a regime where both the temperature T and the size of the system L are finite [38]. 
In particular, one is interested in knowing if the transition between the two very distinct 
behaviors found for the infinite system is still attainable and, if so, how the critical value 
of the pertinent parameter is going to be affected when T and L are finite. To grasp how 
the results obtained so far in this section are going to be changed due to these finite size 
effects we look at the behavior of the two point correlation functions of free fermions on 
the torus [39]. For this system one has that, 



tti 



v\T It) 



(4.5) 



v (0\t) Oi(i\u - v\T \t) ' 

where 6 u (uj\t) are the modular Jacobi theta functions [40], cL#i(0|t) = 5 aj 0i(w|r)| w= o, 
v defines the boundary conditions for tp and r = iLT. For instance, for finite temperature 



- 13 - 



boundary conditions, only the sectors v = 3, 4 of the spin structure of the fermion contribute 
(4.5); this is to say that, on the torus, one can only choose for ip either v = 3, corresponding 
to antiperiodic-periodic (Neveu-Schwarz, NS - Ramond, R) boundary conditions, or v = 4 
which corresponds to antiperiodic-antiperiodic (NS-NS) boundary conditions. 

When LT — > oo, using the standard representation of the 9 U functions [40, 41], one 

gets 

(i/>(u)il>(v)) m = 4ginh J r[M _ ^ [1 ± 2e-" LT cosh27rr|u - v\ + ■■■] , (4.6) 

As a result, in the limit of finite T with L> \u — v\, one may approximate Eq. (4.3) with 
A = 1/2 in terms of (4.5) and write (4.1) as, 



l(A:B) 



c [" T^jui) T(u 2 ,vi) 
3 [T(vi,u 1 )T(v 2 ,u 2 ) 

where T(u,v) is given by [39] 



T(it, v) = log ■ 



- v 2 \T \t) 6 v {i\u 2 - v\\T |t) 
%{i\ui - v\\T \t) 9 v (i\u 2 - v 2 \T \t) 

dM0\r) 



(4.7) 



6\(i\u — v\T \t) 

For L 3> \u — v\ one has that T(u, v) ~ log l/(i|u — v\T); as a result one has 



(A'-B) 



log 



1 



+ |log 



.(il^i - v 2 \T |r) Q v {i\u 2 - v\\T |r) 



(4- 



(4.9) 



,(i\ui - vi\T |r) 6 u (i\u 2 - v 2 \T |r) 

where x is defined in (1.13). One notices that (4.9) reduces to (1.14) when the separation 
between the intervals is very small since the function /„(x,r) defined as 



U(x,t) = log 



,{i\ui - v 2 \T \t) 6 v (i\u 2 - v\\T \t) 



(4.10) 



,{i\u\ - vi\T \t) 9 u (i\u 2 - v 2 \T \r) 

approaches zero when x — > 1 while f u (x,r) > for x < 1. 

We observe that, in a rather large range of values for T and L, there is still a transition 
between two very different behaviors of the MI. However, the critical value xo, at which 
the transition occurs strongly depends on the ratio l/L i.e, l/|r| = (LT)" 1 as reported in 
Fig. 4. Indeed, a numerical analysis shows that, for a finite system, xo is always xq < 1/2 
and that, only as L — > oo, xq — > 1/2 recovering the result in [19]. 

For the sake of completeness we shall compute the quantum correlators and the MI in 
a rotating BTZ black hole background in Appendix A. 



5 Concluding Remarks 

Originally developed within string theory, the AdS/CFT correspondence provides a geo- 
metrical framework to investigate also strongly coupled condensed matter and spin sys- 
tems at criticality . An intriguing observation has been that MERA [25] may be efficiently 
described through the AdS/CFT correspondence by introducing an AdS metric in a perti- 
nently engineered bulk space [22, 24]. In this paper we use this MERA induced AdS/CFT 
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1/|T| 

Figure 4. Dependence of the transition point xo on |r| = LT i.e l/|r| oc l/L for the non rotating 
BTZ black hole (Eq (4.9), circles) and for the quasi-extremal rotating BTZ black hole (Eq (A. 3), 
squares). 

correspondence to provide a framework in which the mutual information and the two point 
quantum correlations between disjoint blocks of a quantum system at criticality may be 
evaluated. We feel that an advantage of this approach is that, at least in principle, is not 
strictly confined to the analysis of one dimensional critical systems. 

In order to get an hint on the pertinent metric to be used to describe the MERA 
induced bulk AdS^+2 space, we observed here that, when computing the quantum corre- 
lations between two disjoint blocks of a boundary quantum critical system, the structure 
of the causal cones characteristic of MERA implies the existence of two different regimes 
attainable by tuning the ratio between the size and the separation of the disjoint blocks. 
To account for this transition we proposed that the MERA induced holographic dual bulk 
spacetime could be described by an AdS^+2 black hole and used the RT formula to compute 
the MI of two disjoint regions of the boundary critical system. Intuitively speaking, this 
amounts to orbifolding the AdS geometry introduced in [22, 24] when dealing with disjoint 
blocks. 

As an explicit example, we used a BTZ AdS3 black hole to compute the MI and the 
quantum correlations between two disjoint intervals of a one dimensional boundary quan- 
tum critical system: here, our analysis not only confirmed the existence of the phase transi- 
tion emerging when the conformal four point ratio reaches a critical value but also provided 
a rather intuitive entropic argument explaining the source of this instability. Furthermore, 
we showed how non universal behaviors may emerge in the holographic computation of the 
MI between two well separated disjoint blocks. Of course, our analysis does not exclude the 
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possibility that other geometries -such as Lifshitz geometries- may account for the behavior 
of the causal cones of disjoint blocks in MERA. 

A remarkable feature of the RT approach to the computation of MI and EE taken in this 
paper is that it associates with each spatial region of the boundary a unique spatial region 
of the bulk [6]. This bulk to boundary map -via the structure of the causal cones- seems 
to play an intriguing role also in MERA. Indeed, we exploited this map in MERA to give 
an ansatz for the dual holographic geometry associated to a region made of two disjoint 
blocks of the (f-dimensional boundary critical system. We, then, observed that -when 
the separation between the two blocks exceeds a critical value- the quantum correlations 
exhibit a thermal behaviour and the EE may be computed as the thermodynamic entropy 
associated to a certain black hole. In the context of the AdS/CFT correspondence thermal 
states have been recently constructed in [7]. 

The AdS/CFT correspondence, is a strong-weak duality. This amounts to say that, 
when the dual gravity description of a quantum system is classical, the correlations on the 
boundary theory are quantum and highly non-local (entanglement) [28]. In the MERA 
induced AdS/CFT correspondence, the locality of the emerging AdS space is due to the 
existence of entanglement at all scales in the quantum critical system located at the bound- 
ary. We feel that our results may help to elucidate the nature (quantum and/or classical) of 
the correlations computed using the RT formula within the AdS/CFT correspondence. In- 
deed, despite the fact that MI quantifies both classical and quantum correlations, recently, 
in [42], it has been shown that MI, when computed using the holographic RT formula, 
obeys the same monogamy relations required for a true measure of entanglement. Since 
the monogamy relations severely limit the amount of entanglement sharable between the 
different parts of an arbitrarily partitioned system [43] this should imply a truly quantum 
nature of the correlations measured in the holographic computation of the MI. In [15], 
using numerical methods to compute a true measure of entanglement such as negativity, it 
has been found that the entanglement between disjoint intervals in spin chains at criticality 
also showed a crossover from pure algebraic decay to pure exponential decay when a critical 
ratio between the separation and the size of the intervals was reached. 

Finally, we feel that the use of a pertinent metric in the AdS space built from the 
MERA induced AdS/CFT correspondence may be exploited also as a way to look for 
alternative and- hopefully- more powerful ways of optimizing MERA tensor networks. 

A Quantum correlators and MI in the rotating BTZ background 

In this appendix we compute the MI and quantum correlations between two disjoint inter- 
vals of the same size I when the background metric is a rotating BTZ black hole. 

When the distance between the two intervals is small enough, the geodesic of minimal 
length is L^Jg = L{u\,v 2 ) + L{v\,u 2 ); using (3.13), one gets 

sinh 2 (itTlqi) sinh 2 (tvTrQi) 



!{A:B) = 3 



sinh (ttT l q 2 ) sinh (ttT l q 3 ) sinh {itT R g 2 ) sinh (ttT r q 3 )_ 

where q\ = \u\ — v\\ = \u 2 — v 2 \ = I, g 2 = \u 2 — v\\ = d and £3 = \u\ — v 2 \ = 21 + d. In 
equation (A.l), the event horizon is located at z + = l/r+ = I and, upon introducing the 
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two variables = l/(r + + r_) and z/j = — r_) such that zi < z + < zr, one is able 

to define the two temperatures Tl = 1/2ttzl and Tr = 1/2ttzr (Fig. 5). 

For the near extremal BTZ black hole, i.e, for a black hole in which Mi > J (r+ > r_, 
zr — >■ oo), equation (A.l) may be written as, 



c , / j; \ c 



/ (^) = 3 1 ° g + 3 l0g 



sinh 2 (ttT l qi) 



sinh (nT L (7rT L g 2 ) sinh (itT l q 3 ) 



(A.2) 



One sees from equation (A.2) that the decoupling of the right and left moving sectors 
induced by the presence of two horizons, plus the near extremality condition of the spinning 
black hole, yields an expression for the MI which decomposes in two terms: one -identical 
to (1.14)-depends only on c and the conformal ratio x- and a second identical to (4.1). It 
is easy to convince oneself that equation (4.9) when T = Tl reproduces the second term of 
(A.2). As a result, the MI of two disjoint intervals in a finite system of total length L may 
be written as, 



T 2c , ( x \ . c , (0y{i\u 1 -v 2 \T L \T)e u {i\u 2 -v 1 \T L \T)\ 

kA - B) = T ^ ) + 3 l0g Wi l \m-v 1 \T L \r)eMu.-v 2 \T L \r)) ' (A " 3) 

with r = iLTi,. 

For the rotating extremal BTZ black hole, numerical simulations show that l/|r| 
weakly affects the location of the transition point xq for the MI (Fig. 4). Furthermore, it 
appears from our results and those presented in [19], that the MI is parametrically small 
at xq. As a result, due to the inequality (1.11), one should expect also here a transition 
for the quantum correlations. 

For the rotating BTZ black hole, the behavior of the quantum correlations between 
points located in different disjoint intervals is given by, 



(0(xa)0{x b )) 



PlPr . , (A X A - x B \\ . , {ir\x A - x B \ x 
— ^— sinh sinh 



h J V fa 



(A.4) 



When the near extremality condition is satisfied, i.e, when r_|_ > r_, so (3r — > oo, the 
two point correlations behave as, 

( (X J)0 ( IS )) J ' V A f^V'expf-M^) . (A .5) 



\x A -x B \J exp \ fi L 

Eq. (A. 5) shows the existence of a crossover from pure algebraic decay to pure expo- 
nential decay of the quantum correlations as fx = \xa — xb\/z+ increases. 
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